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Abstract. We compute the correlation function in the equilibrium version of Renyi’s 
parking problem. The correlation length is found to diverge as 2“^7r“^(l — p)~‘^ when 
p 1 (maximum density) and as 7r“^(2p — 1)“^ when p \ 1/2 (minimum density). 
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1. Introduction 

A parking configuration of hard rods (or cars) on a line or a circle is a configuration where 
the largest gap between successive rods is less than the length of a rod: the parking is 
full. The statistical properties of parking configurations obtained by random sequential 
parking (or adsorption) were first studied by Renyi [1]. The present paper deals with the 
corresponding equilibrium model. 

The correlation function of a gas of hard rods in one dimension was computed by 
Frenkel [2]. A general one-dimensional fluid with a nearest neighbor interaction, strongly 
repulsive at short distance and decaying rapidly at large distance, was then solved by 
Giirsey [3] in a grand canonical ensemble. Salsburg, Zwanzig and Kirkwood [4] derived 
a similar solution in the canonical ensemble, suitable for testing the then newly invented 
Kirkwood-Salsburg equations. 

We give a solution which is a little simpler, in a canonical isobaric ensemble, and 
give precise asymptotic forms for the diverging correlation length near maximum density 
(pressure going to -|-cxd) and near minimum density (pressure going to — cxd). 

The intervals between hard rods can be mapped onto the gradient variables in a one¬ 
dimensional interface model. One could prove mathematically the equivalence of ensembles 
for the parking model, including uniqueness of correlation functions, using the local central 
limit theorem like in the proof of the Wulff shape for one-dimensional interfaces [5]. 

As a by-product, we give the probability for N points distributed uniformly and 
independently in (0, L) to be in a parking configuration, as a function of p = N/L, in the 
thermodynamic limit. 
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2. The model 


A parking configuration of N hard rods (or arcs) of length one on a circle of length L is 
specified by N positions Xi,.. .X^ G [0, L) = R/LZ with 


\x,-x,\>i yi^j 

niin(Xj — Xj)j^ <2 Vi 
niin(Xj — Xi)^ <2 Vi 


( 2 . 1 ) 


The Lebesgue measure on the set of parking configurations, a Borel subset of R^, nor¬ 
malised by the total measure of this set, defines a probability measure This 

probability measure inherits the rotation invariance of the Haar measure on the circle of 
length L. The canonical partition f un ction Zjv,l is defined as 

Zn,l = ^ dxi...dxN ( 2 . 2 ) 

’ • J Parking 

where the range of integration “Parking” is defined by (2.1). The free energy is F{N, L) = 
— IuZn^l- Temperature plays no role and is omitted. 

For x,y ^ [0, L) we define 


N 


PN,Li^) =Y1 

i=l 

N 

PN^,Li^^y)= 

i^j = l 


^N,L{Xi e{x,x + dx)) 


dx 


^N,LiXi e {x,x + dx), Xj e{y,y + dy )) 
dxdy 


(2.3) 


By rotation invariance p^j^\{x) is independent of x. Since 

[ dxp%\{x)=N, 

denoting N/L = p, we get p^^\{x) = p. Conditions (2.1) imply 1/2 < P < 1. Similarly 
Pn\(^^ y) depends only upon \x — y\, and 

[ dx [ dy pP^{x,y) = N{N-1). 

Jo Jo 

The pair distribution function is defined as 


9n,l{x) 



(2.4) 
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and satisfies 


dxgN,L{x) = (l - —) L 

The degree of correlation or independence between the conhgnration aronnd x and the 
conhgnration aronnd y may be estimated by looking at the pair correlation fnnction 



FN,L{^i : XjE {x,x + dx), 3j : Xj e {y,y + dy)) 
Pjv,l( : Xi e {x,x + dx) ) IPw,l( : Xj e {y,x + dy) ) 


gN,L{x -y)-l 


3. S pacings 

Let (Xi,^ 2 ,..., X'j^) be obtained by permntation of {Xi,... ,Xn) so that 

Xi<X'2<...<X'j^<Xi + L 

and let 

Si =X!^- Xi 

S, = X'^,-Xi i = 2,...N-l 
Sn=Xi + L-X'^ 

Then (Xi, Si ,..., Sn-i) is distribnted according to the Lebesgne measnre on the subset 
of dehned by 

0 < Xi < L 
l<Si<2, i = 

N-l (^-1) 

l<L~^Si<2 


and we may write (2.2) as 

Zn,l = j dxi J dsi... J dsjv-i 

= Ly dsi... J dsN 5 Sj — L 


(3.2) 


It follows that (iSi ,..., Sn) is distributed according to 

Fn,l {Si G (si, Si + dsi),..., Sn G (sjv, sn + ds^) ) = 
N N 

n dsi )5f Si — L\ j noi'malisation 


(3.3) 


i=l 
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Then, for a: > 0, 


IPAr,L( j : Xi G (0, dy),Xj e {x,x + dx)) = N'S’n,l{Xi G (0, dy), 3j : Xj G {x,x + dx )) 


N-l 


N ^ ^ Pjv,L (G (0, dy), ^1 + 5*1 + ... + Sm G {x, x dx )) 


771=1 

N-l 


^ PAr,L( Si + . . . + Sjn e {x, X + dx) 


m=l 


so that 


9n,l{x) 


N-l 


E 


^N,l{Si + . . . + Sm ^ {x,X + dx)) 
pdx 


and gN,L{-x) = gN,L{x). 


a: > 0, 


(3.4) 


4. Isobaric ensemble 


In the thermodynamic limit, L —> cx), N ^ oo with p = N/L hxed, asymptotic statistical 
properties are easier to compnte in an ensemble where the pressnre p is hxed instead of 
the system size L. This is dehned as follows: for p G M, let 5'i,...,5'jv be distribnted 
according to 


N 


^N,p{Si G (si. Si + dsi),..., S'at G (siv, sn + dsN )) = ( H dsi j / 

The normalisation is 


norm. 


7=1 


Zn,p — 


ds e 


—ps 


N 


-p -2p 


P 


N 


(4,1) 


(4.2) 


and the associated potential is the Gibbs potential G{N,p) — — In Z^ p. The thermody¬ 
namic relation 

is the nsnal one, with Eiv,p Si = L, the mean system size. We thns get 

^ dG N e~^ - 2e~^^ 

L = — = - \-N -^ 

dp p „ p „ 


or 


e — e 


11 ~P r, -2p 

1 1 e — 2e 


P e ^ - e 


(4.3) 
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Although temperature plays only a dummy role, one may wish to have it in, and also to 
replace the length one of rods and maximum allowed gap between neighboring rods by a 
length 1. Then one gets 


1 

P 


kT 

-h 

P 


—pl/kT —2pilkT 

e — 2e 


^-pljkT _ ^-2pl/kT ’ 


(4.4) 


which looks more like an equation of state. 


1 / 2 . 


Going back to (4.3), for p 


+CX), In 2, 0, —oo, we have respectively p 


1, ln2, 2/3, 


In this ensemble we can compute, for a: > 0, 


N-l 

9N,pi^) = 

m=l 

N-l 

= 5 : 

m=l 

A^-1 

= 5 : 

771=1 

N-l 

= 5 : 

777=1 


^N,p{ Si + . . . + Sm ^ {x,X + dx) ) 


P 


—p —2p 

e — e 


pdx 

m «2 


P 


—p —2p 

e — e 


^ e e dSjnS(y^Si 

j dtl ■ ■ ■ j 5 ^ ti — X + m 


X 


m -px 

e 


P 


—p —2p 

e — e 


m -p^ 

e 


Um{x - m) 


where the densities Un{x) satisfy the recursion relation 


(4.5) 


UnN^) ^ / Un{x-y)dy, 

Jo 


which gives [6, p. 27] 

'Ui(a: — 1) 

Um{x - m) 

The thermodynamic limit g{x) = liniN^oo 9 N,p{x) is given by (4.5) (4.6) with the range of 
m extended to infinity. Note that for each x only [a:/2] terms contribute, because Um{x—m) 
vanishes for x > 2m. Plots of 9 {x) are given in Figure 1. 

IPjv,p(-) also induces a distribution for parking configurations of N rods on [0, cx)), 
so that 9 n,p{x) may be considered to be defined by (2.3)(2.4) with L replaced by p and 
1/p = En,pL/N. Then gN^p{-x) = gN,p{x) and g{-x) = g{x). 


= 1 


xe{i,2) 

1 


m 


im — 1 )! 

^ ’ 1=0 


{-Y 


X 


m 


m— 1 


m > 2 


(4.6) 
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Fig. 1: The pair distribution function g{x) 


5. Correlation length 


Let us compute for A > 0 the Laplace transform, 

-Ar 

^Laplace ^^ dx 

Jo 

-(A+p) -2(A+p)) (5.1) 

p e — e 

~~~P -iX + p){e-P - e-"^) + 

Using (4.3), we find 1 as A —0, in agreement with g{x) —1 as a: —cx). 

We then consider the Laplace transform of the pair correlation function h{x) = g{x) — 1, 


^Laplace^^^ _ ^Laplace^^^ 


1 

A 


and look for its poles, which can only be in the complex half-plane Re{X) < 0. We must 
solve 

-(A+p) -2(A-|-p) -p -2p 


e — e 

X + p 


e 


— e 

p 


( 6 . 2 ) 
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Setting \ = a + ih and p = p + a, the modnlns and the phase in (5.2) give 

/ -P -2px2 . -p -2px2 

(e — e j + 2e (1 —coso) (e — e j 

~ (5.3) 

sin6 —e sin 26 6 

cos 6 — cos 26 P 

When p +CXD, looking only in the strip —1 < a < 0, we get 6 from the second eqnation 
and then a from the hrst: 

b = 2n7r{l-^ + 0{^)) 

99 (5.4) 

2n^7r^ , \ 

where n G Z \ {0}. When p —> — cx), we get similarly 

4 = „.(i_i. + o(l)) 


2 2 
n^TT^ 


a = —- 


(i + o(i)) 


The pair correlation fnnction /i(a:) is conveniently recovered by Fonrier transform: 

7/rr ^ ^ 

e g{x) dx = 

-OO 

I -ikx^, , 


h{x) = 


— ikx ^ 

e h[k) dk 


Using (5.1)(5.4)(5.5), the contonr of integration can be shifted past the nearest poles, 
whose residnes are the dominant part of h{x) when x oo, giving the inverse correlation 
length 

{ 27r^ / 

^ 5 -(l + C>(-)) as p ^ +00 

b 1 


or, nsing (4.3), 


4(2p-l)-" 


as P 1 
as p \ 1 /2 


The imaginary part of the poles gives the psendo-period of oscillation. Ax = 1 when p y' 1 
and Ax = 2 when p \ 1/2, in accordance with Fignre 1 and with the Dirac train limits 



✓ —A 

e 

as 

p - 

■9 +CX) 


lim ^Laplace ^ 

p—i-±CX) 

l-e 

-2A 

2e 

. 1 -2A 

^ 1—e 

as 

p - 

■9 —OO 

(5,9) 
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lim g{x) = I 

p^±oo 2 ^ 5{x-2n) 

^ n e Z 


5{x — n) as p —+cx) 

n G Z 


as p ——cx) 


(5.10) 


6. Equivalence of ensembles 


For any p G M we may write (3.2) as 


Zn,l 

L 


= j dsi ... J dsN Sj ~ L 

= ^Af.P Eaj.p 5 ( ^ — l) 


dJ2st-L 


where Ejvp(-) is defined by (4.1)(4.2). Let ns choose p so that 

IEAr,p ^ Si = L 

where L is the fixed valne in the fixed L ensemble. Then the Central Limit Theorem 
implies 

and we get 

^ Zn,p 

The free energy per particle at density p, 

M = ita 

J N^oo N 

and the Gibbs potential per particle at the corresponding pressnre p, 

,(,)= (6,2) 

N —^oo iV p 


are therefore related by 


f{p) = --+9{p) 
P 


The local version of the central limit theorem can then be nsed to prove 


lim gN,N/p{x) = lim gN,p{^) = 9 {x) 

N^OO N^OO 



along the lines of the proof of the Wnlff shape for one-dimensional interfaces in [5]. Eqniva- 
lence of ensembles is of conrse very standard, bnt its jnstihcation or derivation is generally 
more involved than for the present model as ontlined above. 


The explicit form of the free energy f{p) allows to answer an old qnestion in statistics 
[7,8]: let Xi,..., Xiv be independent random variables each distribnted nniformly over 
the interval (0, L). Denote the corresponding probability distribntion (the ideal gas), 
and call Parking the event that the smallest and largest spacings are respectively larger 
than one and smaller than two. We have 


(Parking) 


dsi... dsjv d(J2si-L) 

Jq dsi... dsN L) 


Zn,l/L 


Let 1/2 < p < 1. Then (6.1)(6.2)(6.3)(6.4) and Stirling’s formnla give 


lim 

N —^oo 


-P _ -2p 

Jj In pjj'y p(Parimg) = 1 - ^ - In p - In ^ - 


(6.4) 


(6,5) 


where p is related to p by (4.3). 
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